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Abstract. Oscillatory instabilities in Hamiltonian anharmonic lattices are known 
to appear through Hamiltonian Hopf bifurcations of certain time-periodic solutions 
of multibreather type. Here, we analyze the basic mechanisms for this scenario by 
considering the simplest possible model system of this kind where they appear: the 
three-site discrete nonlinear Schrodinger model with periodic boundary conditions. 
The stationary solution having equal amplitude and opposite phases on two sites and 
zero amplitude on the third is known to be unstable for an interval of intermediate 
amplitudes. We numerically analyze the nature of the two bifurcations leading to this 
instability and find them to be of two different types. Close to the lower-amplitude 
threshold stable two-frequency quasiperiodic solutions exist surrounding the unstable 
stationary solution, and the dynamics remains trapped around the latter so that in 
particular the amplitude of the originally unexcited site remains small. By contrast, 
close to the higher-amplitude threshold all two-frequency quasiperiodic solutions are 
detached from the unstable stationary solution, and the resulting dynamics is of 
'population-inversion' type involving also the originally unexcited site. 
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1. Introduction 

As indicated by a number of recent results (e.g. [HEllSllllEllSllZllSlin]), oscillatory 
instabilities of time-periodic solutions in Hamiltonian anharmonic lattice models are 
quite ubiquitous, and likely to play an important role in the dynamics of any physical 
system which, in some approximation, can be described by such a model. Typically they 
are associated with an inhomogeneous amplitude distribution of the solution dividing 
the lattice into two sublattices of small- and large-amplitude sites, respectively, and 
oscillatory instabilities appear when internal oscillations of the sublattices resonate. 

As long as the harmonic coupling is not too strong, time-periodic solutions for 
coupled anharmonic oscillators can be conveniently described using the multibreather 
formalism ^Hl utilizing the ('anticontinuous') limit of zero coupling to identify 
each solution with a coding sequence. In the simplest case of time-reversible solutions 
with one single determination of the oscillation pattern for the individual uncoupled 
oscillators at a given period, the code is defined as +1 for each site oscillating with a 
certain reference phase, -1 for each site oscillating with the opposite phase, and for 
each site at rest. (This scheme can be readily generalized to include e.g. higher-harmonic 
oscillations by allowing codes ±2, ±3, or non-time-reversibility by allowing arbitrary 
phase torsions [HI). This identification is unique for sufficiently weak coupling [TUl[TT| . 
and can be obtained by numerical continuation using Newton-type schemes (e.g. [T^). 

In particular, for Klein-Gordon (KG) chains of harmonically coupled soft 
anharmonic oscillators, most known examples of oscillatory instabilities are associated 
with solutions whose coding sequences contain blocks of the kind 1, 0^, —1, 0, with 
integer p > (0^ means repeated p times). (For hard oscillators, a corresponding 
statement is true provided that the code sequence is multiplied with (—1)".) Thus, 
large-amplitude sites of alternating phases (codes ±1) and small-amplitude sites (code 
0) provide a two-sublattice structure. Generically, such solutions are linearly stable 
for weak coupling when the internal oscillation frequencies for the two sublattices are 
separated, but for larger coupling the frequencies may overlap causing resonances and 
instabilities. By contrast, blocks of the kind 1, 0^, 1, 0, ... typically yield nonoscillatory 
instabilities already for weak couplings. This difference can be related to the different 
natures of extrema for a corresponding effective action, defined in jTT] as a function of 
the individual phases at each site (see also [121 [1110], and [T3] for further developments). 

Oscillatory instabilities were explicitly demonstrated e.g. for: (i) two-site out-of- 
phase breathers |I| [Ilj ('twisted localized modes' [ISIISl) with code ...0,0,1,-1,0,0...; 
(n) discrete dark solitons [HI HZI [Zj ('dark breathers' [Hj) with a hole and a phase shift 
added to a stable background wave, code ...1,— 1,1,— 1,0, 1,-1, 1,-1...; (in) nonlinear 
standing waves [11[Z1[H1, corresponding to periodic or quasiperiodic repetitions of blocks 
of codes of type {1, 0^, —1, 0''} with particular combinations of integers p, g > 0. These 
instabilities may play an important role in the processes of breather formation and 
evolution to energy equipartition [Zl[Hl- Oscillatory instabilities are found also e.g. for 
two-dimensional vortex-breathers ^3 ^1 [0] . 
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Mathematically, oscillatory instabilities appear through 'Krein collisions' (e.g. [TU] ) 
of eigenvalues to the corresponding linearized equations of opposite Krein (symplectic) 
signatures^. Although well-known from other branches of physics, their relevance for the 
dynamics of anharmonic lattices was realized mainly after |TI] . Generally linear stability 
of time-periodic solutions is studied through the linearization of a Poincare return map 
with the exact solution as a fixed point (Floquet map), and its symplectic nature implies 
that all eigenvalues (Floquet multipliers) must lie on the unit circle for stability. In a 
Krein collision, four eigenvalues leave the unit circle as a complex quadruplet. 

However, for weakly coupled KG lattices the small-amplitude dynamics can be 
approximated, through a multi-scale expansion (e.g. [211 IZ|), by a discrete nonlinear 
Schrodinger (DNLS) equation. A time-periodic KG solution then becomes a purely 
harmonic (stationary) DNLS solution, which can be made time-independent by studying 
the dynamics in a frame rotating with this frequency. Thus, performing also the stability 
analysis in the rotating frame reduces the problem to studying stability of (relative) 
equilibria, as no explicit time-dependence appears in the linearized equations |221I2S1- 
Linear stability demands all eigenvalues of the linearized flow to lie on the imaginary 
axis, and a Krein collision manifests itself as two pairs of complex conjugated imaginary 
eigenvalues with opposite Krein signature colliding and going out in the complex plane as 
a quadruplet (e.g. |211)- This collision is also (e.g. [2312^1) referred to as a 1:-1 resonance 
(equal frequencies, opposite Krein signatures), and corresponds to a bifurcation of two 
(relative) time-periodic solutions which, due to its similarity to a Hopf bifurcation in 
dissipative systems, has been termed [2Z] a Hamiltonian Hopf (HH) bifurcation. 

Following a recent description of the essential aspects of the HH bifurcation ^| , we 
briefly recall its phenomenology (see 123121] for a rigorous description). HH bifurcations 
may be of two different types called [2H1 'type L and 'type H', respectively (in a 
related context, these bifurcation classes were called [221 'hyperbolic regime' and 'elliptic 
regime', respectively, referring to the nature of the corresponding normal forms). In the 
type I bifurcation, two ('Lyapunov') families of stable periodic solutions surround the 
equilibrium in its stable regime, merge at the bifurcation point, and form a single family 
detached from the equilibrium in its unstable regime. For type II, the two periodic 
solutions belong to a single connected family on the stable side of the bifurcation, which 
shrinks to a fixed point and disappears on the unstable side. Thus, while for type I stable 
periodic solutions exist close to (but away from) the equilibrium also in its unstable 
regime, for type II no periodic solutions exist close to the unstable equilibrium. Since for 
the DNLS equation this description refers to the rotating frame, the (relative) equilibria 
correspond to time-periodic stationary solutions and the (relative) periodic solutions 
generally to quasiperiodic two-frequency solutions in the original frame. For the DNLS 
equation (and other equations with particular symmetries [SD]) such quasiperiodic 
solutions may persist, and even be spatially localized, also for infinite lattices [THl El [Sl- 
it is our purpose to describe, mainly by numerical means, the nature of these 

X We here conform to standard terminology, although historically the concept of symplectic signature 
was known before Krein, see e.g. [201 • 
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bifurcations, aiming at an increased understanding for the dynamics in the oscillatory 
instability regime. We focus here on the simplest possible model where the above 
described mechanism for oscillatory instabilities can occur: the 3-site DNLS-model with 
periodic boundary conditions and the stationary solution with code +1, —1, 0. The fact 
that this solution exhibits a complex instability is known already since the pioneering 
works of Eilbeck et.al. [221123; however, in spite of the wealth of existing literature on 
this model (e.g. |221l23EIlElESllSlESlEniEI| and references therein), apparently 
the unstable dynamics of these stationary states and the properties of the surrounding 
quasiperiodic solutions were never before analyzed. Future work will extend our results 
for the DNLS trimer to the above-mentioned examples for infinite lattices. 

In addition to its importance for general anharmonic lattice dynamics, the DNLS 
model has many other particular applications (see e.g. jHHl IHn] for recent reviews). 
In particular, the recent experimental developments in the fields of coupled optical 
waveguide arrays (e.g. jlOl HI]) and Bose-Einstein condensates (e.g. |121 EHl EZl 1^ ) 
provide very exciting possibilities for practical utilization of the theoretical results. 

The outline is as follow. In section |21 we define the DNLS model and describe the 
stationary {1,-1,0} solution and its linear stability properties, partly recapitulating 
material from j22ll2SI- In section El we use the technique developed in [T31 to calculate 
numerically families of exact quasiperiodic two-frequency solutions surrounding the 
stationary solution, and analyze the nature of the two HH bifurcations associated with 
its Krein collisions. We find them to be of different types, with a type I bifurcation 
at the low-amplitude instability threshold for the stationary solution and a type II 
bifurcation at the high- amplitude instability limit. In section 0] we show that, as 
a consequence, the instability-induced dynamics of the {1,-1,0} solution exhibits a 
new type of self-trapping transition: in the low-amplitude part of the unstable regime 
the dynamics remains trapped around the initial solution, while in the high-amplitude 
part the dynamics also to a significant extent involves the initially unexcited site in a 
'population-inversion' jH] type dynamics. Section El concludes the paper. 

While this work was in progress, we became aware of a preprint of (published 
after our original submission) which, within the framework of describing the dynamics 
of three coupled Bose-Einstein condensates, reports numerical simulations exhibiting 
self-trapping effects and 'population-inversion' dynamics similarly to our observations. 
Whenever appropriate we will discuss the relation between our results and those of [Hj . 

2. The DNLS model, stationary solution and linear stability 

We consider the following form of the DNLS equation 



(1) 



where we assume C > 0. It is the equation of motion for the Hamiltonian 



i/({iV'4,{C}) = 5^ C|^„,+i-^. 



n=l 



n 



(2) 
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with canonical conjugated variables {i^„}, {'>Pn}- The DNLS equation has also a second 
conserved quantity, the norm§ (excitation number) H = J2n=i iV'nP- We consider the 
trimer, = 3, and impose periodic boundary conditions i^o = ipz, ipi = i'l describing a 
triangular configuration with identical coupling between all sites. Note that we consider 
the case of positive intersite coupling and negative (attractive) nonhnearity in Q, but 
the general case can be recovered since changing the sign of C in is equivalent 
to making the transformation ipn i—^)^Q~^'^^^'ipn, while the same transformation 
followed by a time reversal t ^ —t is equivalent to changing the sign of the nonlinearity. 
However, as N here is odd, these transformations must be accompanied by a change of 
boundary conditions from periodic to antiperiodic. Thus, our model is not equivalent 
to the periodic three-site configuration with positive coupling and positive (repulsive) 
nonlinearity considered in jH]. In the latter case, the obtained solutions are equivalent 
to a subclass (antisymmetric under transformation n — >■ n + 3) of the solutions to 
with C > for a hexagonal geometry (A^ = 6) analyzed in We also note that when 
C 7^ 0, it can be normalized to C = 1 by the transformation ijjn — \/Cil)n-i t t/C. 
A time-periodic, stationary solution to (Q) with frequency A has the form 

^i^Ht) = 0i^^e^^* (3) 
with time-independent ^i^^ We consider the particular family of solutions 

{00^)} = {A, -A, 0}, with A = VA + 3C ; A > -3C, (4) 
for which the norm and Hamiltonian can be expressed as 

Af = 2A^ = 2{A + 3C), (5) 

H = A^{6C - A^) =9C^-A^ = -^(qC--^]. (6) 



2 V 2 ^ 

Linearizing (P) around a stationary solution by writing ipnit) = (^<pn^^ + en{t)j e'^* yields: 

ie„ + C(e„+i + e„_i - 2e„) + 2|0(,^f e„ + (0i^))'e: - Ae„ = 0. (7) 

For real (pi!^^ a linear eigenvalue problem is obtained by writing e„ = a„e~"^'* + 6*e"^**: 

(2C + A - 2(0(f))2)a„ - C(a„+i + a„_i) - {<Pl,^^)%n = ^ia„ 

(0i^))2a„ + C(6„+i + 6„_i) -{2C + A- 2(0i^))2)6„, = oo^b^. (8) 

Then, linear stability is equivalent to all eigenfrequencies uJi being real. Alternatively, 
one may use the linear combinations f/„ = a„ + Wn = a„ — 6„ so that ^ becomes 

jCoW^ = {2C + A- (0(f ))')W^„ - C{Wn+l + Wn-l) = UJ,Un 

£it/„ = (2C + A-3(0i^))2)f/„,-C(f/„+i + [/„_i) =uj,W^. (9) 

§ We here conform to established nomenclature in the physics literature, although mathematically Af 
is rather the square of a norm. 
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The usefulness of this formulation becomes clear by noting that, if e„ is expanded in its 
real and imaginary parts e„ = e"n + ien"*, (0) can be written in the Hamiltonian form 

{4"} )-\-c, <>){ {.<?} j <^°) 

where £o and £i defined by Q are symmetric. Then, if iuj\ is a purely imaginary 

eigenvalue of (|l(Jj) with c^i > 0, the corresponding eigenvector is [ ) where 

{Un, Wn) is the (real) solution of Q with eigenfrequency ui. Thus, to each real 
eigenfrequency c<Ji > of (jSj) and Q we can associate a Krein signature /C(co'i) from 
the symplectic product of the corresponding eigenvector with eigenvalue icji of (fTUI) with 
itself, which, using the standard sign convention as in e.g. reads 

/C(^i) = signJ^UnWn = signJ2 ' M ■ (H) 

n n 

The Krein signature is then the sign of the (Hamiltonian) energy carried by the linear 
eigenmode [201 • 0-^ earlier publication ^E], an opposite sign convention was used.) 

Since is invariant under global phase rotations, cJi = is always an eigenvalue of 
algebraic multiplicity 2 and geometrical multiplicity 1, corresponding to f/„ = 0, Wn = 

('phase mode'). For the solution (0)), the remaining eigenfrequencies are j22] 



CJi\2 1 /A V 1 



\ 




A\^ , A 



+^c'^^^- ^^^^ 



These are real for small or large A/C, but complex for < A/C < 6 where x^, ~ 1.5385 
is the only real solution to x^ + 6x'^ +4x — 24 = ^[j^ (figure [T]) Close to the linear 
hmit A/C ^ -3 we can expand (HD) using A = A^/C = A/C + 3 = Af/2C as a small 
parameter, and obtain \uJi+\/C ^3 + \/3 + AX^/27 and \uJi-\/C ^ 2A^/V3- For uji+ the 
corresponding eigenmode can be expressed as 

-2A/9 
1 + 2A/3 

with /C(co'i+) = —1, while for the eigenmode is most conveniently described as 





0{\^) (13) 





+ OiX'^^) (14) 

with /C(ci;i_) = +1. On the other hand, close to the anticontinuous limit A/C —>■ 00, 
we obtain from (fT^ the eigenfrequencies ^ A/C + 2 and \uJi^\/C ~ 2^ A/C + 

h^CjA. Then, we can write the corresponding eigenmode for cui^ as 





0{{C/Af) (15) 
with /C(co'i+) = +1, and the eigenmode corresponding to uj\^ as 

+ 0{{C/Af'^) (16) 
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with lC{uj\S) = —1. Thus, the Krein signatures of the frequencies uij^ and uj\- are 
interchanged by the two Krein colhsions at the boundaries of the complex regime, where 
the stationary solution is unstable (and the Krein signature undefined). Note also that 
the eigenmodes corresponding to nonzero Ui are symmetric around n = 3, and thus their 
excitation breaks the spatial antisymmetry of the stationary solution around this site. 

By considering the solution (jU as periodically repeated in an infinite lattice, 
it becomes, in the terminology of [H El EI, a nonlinear standing wave with wave 
vector Q = 27r/3 of 'type H', for which (fT^ gives the subset of all eigenfrequencies 
corresponding to eigenmodes with spatial period 3. Close to the linear limit, (fT^ can 
then be interpreted as a phonon mode with wave vector q = 0, while ()14j) represents a 
translational mode of wave vector q = 2tt/3 corresponding to a sliding towards the 'type 
E' Q = 27r/3 standing wave having a periodic repetition of the codes +1,-1,-1. Close to 
the anticontinuous limit, (fTH|l represents a 'hole mode' localized at the zero-amplitude 
site, while fll6p corresponds to an internal oscillation at the non-zero amplitude sites. In 
the intermediate, unstable, regime, these characters of the eigenmodes become mixed. 

3. Two-frequency solutions: bifurcations and linear stability 

As shown in [T3], exact, generally quasiperiodic, two- frequency solutions to the DNLS 
equation exist and can be obtained from the ansatz 

Mt) = <Pn{ty^'\ (17) 
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so that in the variables 0„ the DNLS equation (Q) becomes 

i0„ + C(0„+i + 0„_i - 20„) + |0„p0„ - u;o0n = 0. (18) 

Thus, the frequency cuq, corresponding to a pure phase rotation, becomes a parameter in 
()18|). and for a given cjq we may search for time-periodic solutions to (fTHj) with frequency 
ciJb fulfilling (pnit + 27r/co'b) = 4>n(t). In particular, any stationary solution to ((T)) of 
the form ^ is trivially a time-periodic solution to p8|l of the form (j)n{t) = 0n^''e"^'"* 
with frequency tUb = A — cuq- However, when the linear eigenfrequencies ui (fT^ are 
real and nonzero, in general (Lyapunov) families of genuine two-frequency solutions 
(generally quasiperiodic in the original frame) bifurcate from the stationary solutions 
when ci;i(A) +ci;o = A. In the limit of small-amplitude oscillations around the stationary 
solution, these have frequencies o^b = wi in the frame rotating with the stationary 
frequency A, and coincide with the linear eigenmodes. 

Using standard numerical continuation techniques described in [Hj, these two- 
frequency solutions can be continued versus c<Jb towards larger oscillation amplitudes 
away from the stationary solutions, and explicitly calculated to computer precision 
in their whole domain of existence. Using standard fortran routines for numerical 
integration and matrix inversion, when necessary in quadruple precision, we may 
generally obtain solutions fulfilling ||0„(t-|-27r/ci;b) — 0n(^)|| < e with e ~ 10~^°, although 
for strongly unstable solutions inaccuracies in the numerical integration might limit the 
attainable precision to larger e. All results presented here have been obtained by first 
demanding a precision e ~ 10^^^, and then checked to be invariant when decreasing e. 

The two-frequency solutions to belong generally to two-parameter families, 
and there is some freedom in choosing independent parameters for the numerical 
continuation. (Since C only yields a rescaling, we put C = 1 in all numerical 
calculations.) One may choose e.g. frequencies wc^Jb, or norm Af and action ('area') 
integral / of the time-periodic solution to (fTHj) . i.e. in the frame rotating with frequency 
uq. Using in this frame the canonical variables {i0n}, {(pn}, corresponding to a 
transformation of the original Hamiltonian into H' = H + uJoAf, I becomes 

n=l ^ n=l 

For a stationary solution with i<pn = —ujh4'n, I is simply proportional to the norm 
(/ = —A/") II, while for genuine two-frequency solutions / and M are independent. In the 
following, we choose as independent parameters the frequency tUb in the rotating frame 
and the norm A/", and perform the continuation versus Wb for fixed norm by, at each 
value of cUb, vary the parameter ujq in (fTH|) so to obtain a solution with the desired value 
of M (this procedure may fail at particular points if the dependence M[ujq) for solutions 
at fixed a;b is nonmonotonous or multivalued; examples of this can be found below). 

II The minus-sign arises only because in terms of action/angle variables, frequencies are defined with a 
sign opposite to that of l|2Jl and ((T7|l . see e.g. ,14, and section^ 
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X 




(21) 



Linear stability of two-frequency solutions is obtained by linearizing (jl8|) around 
an exact time-periodic 0„(t) = 4'n\t) + i0n''(t) {4>n\ 4>n real). The linearized equation 
for the small perturbation e„(t) will then have the same form as ((Zj) with A replaced by 
ciJo, and the real and constant replaced by a generally complex and time-periodic 
<f)n{t)- Decomposing e„(t) = en\t) + ien\t) as in (fTn|) . the linearized equations can be 
written in Hamiltonian form as 

coo - 3(0^:^)2 - 1 - CA -20^02^1 \ / {e«} 

c^o - (0^:^)' - 3(0«)2) 1-CA )[ {e^^} 

(20) 

where 1 denotes the A^-dimensional unit matrix and A the discrete Laplacian, Ae„ = 
e„+i + e„_i — 2e„ (note that the second matrix in (j^Uj) is symmetric). Integrating (j^Ull 
numerically over the period Tb = 27r/u;b of 0^ yields the symplectic 2 x 2 A^-dimensional 
Floquet matrix Fti, defined by 

{el:\T,)}\^ ( {e«(0)} 
{ei'^Tb)} ) [ {e«(0)} 

and the solution 0„(t) is linearly stable if and only if all eigenvalues a of F^^^^ lie on the 
unit circle. In the particular case of a stationary solution, the eigenvalues a of Ftj^ are 
related to the linear eigenfrequencies cui of (jH)), Q as a = e'^'^'^'/'^'^ with c^b = A — cuq. 

For a genuine two-frequency solution there are two arbitrary phases, one 
corresponding to global phase rotations which is put to zero at t = in the ansatz 
()17|) . and another corresponding to the phase of the time-periodic solution in the 
rotating frame which we fix by imposing time reversibility, 0*(— t) = 0n(t) (i.e. 
0n''(O) = 0). To each of these phase degeneracies there is a corresponding solution 
to the linearized equations with period Tb, and thus a pair of eigenvalues a = +1 
of Ft^ j^. The solution corresponding to an infinitesimal global phase rotation is 
e„(t) = i0„(t) (i.e. ^{e^[^}, {cn^jj = (^{-0n^}, {0n^}j), while the solution corresponding 
to an infinitesimal translation of the phase in the rotating frame is e„(t) = 0„(t) (i.e. 
^{e^[^}, {ei'''} j = ^{0n''}, {0ri''} j )• (Note that for stationary solutions, these become 
linearly dependent and there is only one pair of eigenvalues at +1.) Thus, for the 
particular case A^ = 3 there is for a genuine two-frequency solution at most one pair of 
Floquet eigenvalues a not identical to +1, so the only way such a solution can change 
its stability is through collisions of the members of this pair either at +1 or at -1. 

Furthermore, due to the invariance of the DNLS equation and its linearization under 
the transformation 0„ — 0^ and the antisymmetry (02'^'' = — 0^^^;03'^^ = 0) of the 
particular stationary solutions, the families of two-frequency solutions bifurcating from 
them will have the additional space-time symmetry 0i(t-|-Tb/2) = 02 (t), 03()f: + Tb/2) = 
03 (t). This implies that the time-periodic solution 0„ with 0i(O) 7^ 02(0) contains 
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only odd harmonics of its fundamental frequency cUb- As a conseqeunce, collisions 
of Floquet eigenvalues at -1 corresponding to period-doubling bifurcations will not 
yield instabilities, which can be seen by decomposing the Floquet matrix as Fx^ = 
Pi^2FTb/2Piw2FTb/2, where is the matrix obtained by integrating over half 

a period and Pi«_,2 interchanges sites 1 and 2. Thus, only collisions of eigenvalues j3 to 
Pi^2FTb/2 at /3 = ±1 can yield instabilities, and they must all correspond to eigenvalue 
collisions at a = +1 for the total Floquet matrix Y^y^ since a = Below, we report 
the numerical stability properties in terms of eigenvalues j3 rather than a. 

3.1. The low- amplitude 'Type I' bifurcation 

3.1.1. Small M/C. We first discuss some properties of the two families of two- 
frequency solutions bifurcating from the stationary solution for A/C < x^,, i.e. M/C < 
9.077..., corresponding to the linear modes with asymptotic oscillation patterns (fT^ and 
dlH), respectively, for small Af/C. Their continuation versus cUb for fixed Af/C = 5 is 
illustrated in figure El For both solutions, the continuation from the stationary solution 
is monotonous in the direction of decreasing tUb, and they are linearly stable for a 
rather large interval of frequencies < ^^i±, but become strongly unstable at small Wb 
(figure 121 (c)). For the (grey) solution corresponding to uji+, the instability is caused 
by a collision of eigenvalues f3 at -1 producing a rapid and monotonous growth of \f3\ 
for small cjb- This behaviour is typical for M/C > 4.5. For M/C < 4.5, (3 attains a 
minimum value and then returns towards the unit circle, and finally through a collision 
at /5 = — 1 the grey solution bifurcates into the homogeneous stationary solution with 
01 = 02 = 03 = — a/A/'/S, which is linearly stable in this regime |22l ISHI- (The 
continuation versus Uf, is typically non- monotonous close to this bifurcation point.) 
At M/C = 4.5 the homogeneous solution becomes unstable [221; analyzing 
the dynamics of the grey solution in the limit cUb — > for M/C > 4.5, we find that 
it asymptotically approaches a (strongly unstable) solution consisting of homoclinic 
connections of unstable and stable manifolds of the homogeneous solution (as further 
discussed in section 13.1.31 see also figure 01 (b) below). 

The scenario for the (black) solution corresponding to cui- is more complicated: 
it first becomes unstable through an eigenvalue collision at (3 = +1, associated with 
a change of nature of the solution as seen in figure |21 (a). This is related to a third- 
harmonic resonance with the grey solution, and the instability appears at Uh close to 
(but not exactly at) cui^/S. The instability first grows when decreasing cUb, but attains 
a maximum and finally, through a second collision at +1, the black solution becomes 
identical to the (stable) grey solution at frequency 3a;b, which becomes unstable slightly 
after. As mentioned above, for symmetry reasons higher-harmonic resonances can occur 
only for odd harmonics. Since ui- decreases towards zero when M/C decreases while uij^ 
remains finite, for smaller M/C the resonances are of higher order. Thus, for example 
when M/ C approaches from above a critical value ~ 3.53 (where uj\+ = 3u;i_) the value of 
cUb where the black solution becomes unstable approaches the linear frequency tUb = cui-. 



Hamiltonian Hopf bifurcations in the discrete nonlinear Schrddinger trimer 



11 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 

Figure 2. Continuation versus lo-^ at fixed norm A/" = 5 (C = 1) of the two- frequency 
solutions bifurcating from the stationary solution Q with A — —0.5 at ujh = wi± given 
by (|12() : black [grey] curves correspond to lo\- ['^\+\- (a) Initial amplitudes 0„(O) for 
n — 1 (upper), n — 2 (lower) and n = 3 (middle curves); (b) Hamiltonian H (jJl; 
(c) real parts of stability eigenvalues /3; (d) frequency ujq of the rotating frame. The 
horizontal lines in (a) and (b) and the dashed line in (d) correspond to the stationary 
solution according to and The large negative eigenvalues not shown in (c) for 
small ojb will continue to decrease rapidly (e.g. at ujh = 0.2, /3 « —10^ for the grey 
solution). For cuf, < 0.90, the black solution is identical to the grey at frequency 3wf,. 



and for smaller Af/ C the third-harmonic resonance disappears and the solution is stable 
until a similar fifth-harmonic resonance appears for smaller uj\^. It is interesting to 
note, that if the numerical continuation is performed 'carelessly', i.e. with larger steps, 
one might 'jump' the regime yielding a pth harmonic resonance and catch instead a 
solution in apparent (but non-smooth!) continuation of the original solution, as if the 
resonance had not existed. This solution can then be continued towards smaller uo^, 
until a (p + 2)th harmonic resonance is encountered, etc. Thus, a cascade of higher- 
order resonances will result, and the scenario is analogous to that previously found for 
higher-order resonances in Klein-Gordon chains in [7j (e.g. figures 12 and 16 in this 
paper) and |1H]. The behaviour shown in figure |21 is typical for 3.6 < M/C < 5.5, 
while for M/C ~ 5.6 two additional eigenvalue collisions at /3 = +1, producing a new 
pair of small stability/instability regimes, appear slightly before the transformation of 
the black solution into the grey Scub-solution. A further increase of N /C yields one 
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more pair of small regimes of stability /instability through a collision at /? = — 1, and for 
M /C > 6.8 the continuation versus cUb becomes nonmonotonous of 'Z-type', where three 
different solutions of the family have the same frequency cUb, in some regime before the 
final bifurcation with the 3ct;ft-solution (similar behaviour was also observed in [7| I48j ) . 

Note also (figure El (b)) that the Hamiltonian of the black solution changes slope so 
that dH/duj)^ = at the first eigenvalue collision at /5 = +1 (but not at the second where 
the solution transforms to the grey third-harmonic solution). Moreover, the additional 
eigenvalue collisions at /? = +1 for M jC ~ 5.6 also produce a new pair of local min/max 
for H{uJi,) at the collision points. Generally (see also several similar examples below), it 
appears that dH/ dujh = for two-frequency solutions at fixed norm could be a sufficient 
(but not necessary) condition for a change in stability through an eigenvalue collision at 

analogous to the wellknown criteria dA/" /dA = for stationary solutions [Uj and the 
more general dX/ dti; = for periodic solutions of Hamiltonian systems with frequency u 
and action I [21]. Indeed, this could be expected since a turning point of H for relative 
periodic solutions at fixed J\f should correspond to a saddle-centre bifurcation. 

3.1.2. Around the bifurcation point. Focusing on the regime of c^b close to uj\±, we 
compare the solutions on both sides of the bifurcation point M /C ~ 9.077 in figure 
01 For 8.95 < N /C < 9.077, the continuation of the grey solution from the stationary 
solution is no longer monotonous versus tUb, but goes first in the direction of increasing 
uj\y from uj\j^ up to a maximum value, and then turns and continues towards smaller 
tUb similarly as in figure |21 As mentioned above, this multivaluedness causes numerical 
problems for the method of continuing for fixed norm by varying c^o at one point (close 
to but different from the turning point), where dA/" /dtuo for fixed cjb becomes infinite. 

At the bifurcation point N/C ~ 9.077 the gap separating the black and grey 
families of solutions in figure El (a), (c), (e) closes (since uj\- = at this point), and 
for larger Af/C as in figure El (b), (d), (f) they constitute two connected branches of 
a single family of two-frequency solutions, now separated from the stationary solution. 
Note that even though the two-frequency solutions are detached from the stationary 
solution on the superthreshold side of the bifurcation (where the latter is unstable and 
uj\± are complex), the distance between the two-frequency family and the stationary 
solution (which can be expressed in some suitable norm) decreases continuously to zero 
when approaching the bifurcation point from above. This is the characteristic property 
of a type I HH bifurcation |2H] , and the scenario in figure El is equivalent to that of 
figure 2 in j2H| (see also figure 2(b) in [2^1 )• The importance of this property for 
the dynamics of the unstable stationary solution will be illustrated in section |31 

3.1.3. Larger M/C. Increasing M/C on the superthreshold side of the type I 
bifurcation, the continued solutions persist as two connected branches of a single family 
detached from the stationary solution (jH), continuously increasing the distance to it. 
The typical behaviour for larger M/C is illustrated in figure 01 Compared to the earlier 
discussed examples we note two main differences. First, there is a small bubble of 
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Figure 3. Same as figure|21(a)-(c) but for TV = 9.07 (left figures) and TV = 9.08 (right 
figures). In (a), (c), (e) the meaning of the black and grey curves is the same as in 
figure m while in (b), (d), (f) they illustrate the two difi'erent but connected branches 
of coexisting two-frequency solutions continued from the solutions in the left figures. 

instability through a colhsion at /3 = — 1 for the black solution in figure 0] (c) for 
cUb ~ 2.9 — 3.0. In fact, such an instability occurs in many regimes also for smaller 
M jC (e.g. for M jC > 5), but it is then very weak and practically invisible numerically, 
while it grows considerably stronger for larger M jC . Second, for M /C > 9.45 the black 
solution no longer transforms into the third-harmonic of the grey solution for small tUb, 
but approaches as c^b ^ a (strongly unstable) solution with the same values of 0n(O), H 
and ojq as the asymptotic grey solution (figure|ll(a), (b), (d)), but with {3 +oo (figure 
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Figure 4. Continuation versus lo^ at fixed norm A/" = 15 of the two-frequency 
solutions continued from those of figures The meaning of the black and grey parts 
is the same as in figure |31(b), (d), (f); other quantities have the same meaning as in 
figure |2] The eigenvalues not shown in (c) for small Wb continue to decrease [increase] 
monotonously for the grey [black] solution; (3 « —10^^ [(3 « 10^^] at cjb ~ 0.5. 



111(c)). Thus, these asymptotic solutions are not identical, although they both consist of 
homoclinic connections of the unstable homogeneous stationary solution 0i = 02 = 03 
(which has a doubly degenerate unstable eigenvalue [Sj so the corresponding stable and 
unstable manifolds can be connected in different ways). Analysis of the dynamics yields 
that the grey solution in half a period contains two stable and two unstable manifolds 
from the homogeneous solution, and they are all associated with two-site localization 
with a third site of small [0^1^ as for t = in figure IH (a) (see figure IHl (b) below). The 
black solution on the other hand contains three stable and three unstable manifolds in 
one half-period, and of these two pairs are associated with one-site localization with two 
sites of smaller |0„p, and only one pair with two-site localization (see figure IHl (a)). 

We also note, by comparing figures El (d) and El (b), that while close to the 
bifurcation point the Hamiltonian for the unstable stationary solution ^ cuts the 
two- frequency solution family on the black branch (figure El (d)), the corresponding 
intersection point appears for larger A/'/C on the grey branch (figure 131(b)). This change 
of branch occurs at A/'/C ^ 10.6, which, as will be discussed in section HI below, is the 
value where the self-trapping transition of the unstable stationary solution is observed. 
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Figure 5. Same as figure |21 but for the two-frequency solutions bifurcating from 
the stationary solution with A = 7 at A/" = 20; grey [black] curves corresponding to 
solutions bifurcating at Wb ~ ^i- [^b = "^i+l- The eigenvalues not shown in (c) for 
small Wb increase monotonously for both solutions {(3 ^ 10® at LO\y = 0.8). 



3.2. The high- amplitude 'Type IT bifurcation 

Turning our attention to the two-frequency solutions associated with the high-amphtude 
Krein coUision for the stationary solution at A/C = 6, i.e. M jC = 18, we first describe 
some properties of these solutions (which are not identical to those involved in the 
low-amplitude bifurcation) on the stable side, i.e. for Af/C > 18. 

3.2.1. Large M jC . The continuation of the two-frequency solutions corresponding to 
linear modes with asymptotic oscillation patterns (fT3j) and (fTBj) . respectively, for large 
AA/C, is shown in figure El for M/C = 20. For both solutions the continuation from 
the linear modes is monotonous in the direction of decreasing uy,, and they are linearly 
stable in some regime uy, < uji± but become strongly unstable for smaller Ub (figure El 
(c)). The stability regime for the black solution is rather small, and it is unstable in the 
whole range where also the grey solution exists. Also here (figured (b)) dH/duh = 
at the instability thresholds occurring through collisions at /5 = +1 (there is also a tiny 
instability for the grey solution caused by collisions at /5 = — 1 for ~ 4.70 — 4.73). 
As cUb — > also these solutions approach homoclinic connections of stable and unstable 
manifolds of the homogeneous stationary solution 0i = 02 = 03, which however are 
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Figure 6. Continuation of tlie solutions in figure[51for M = 18.5. The behaviour for 
uJb < 5.5 is qualitatively the same as in figure|Sl 

distinct both from each other and from those described in section 13.11 They both 
contain one stable and one unstable manifold in each half-period, but while for the 
black solution the dynamics is of one-site-localization type with two sites of small equal 
|0„P (figureini(c) below), the dynamics of the grey solution have all three sites of unequal 
amplitude as seen in figure El (a) (see figure IHl (d)). 

3.2.2. Around the bifurcation point. Approaching from above the bifurcation point 
N /C = 18 the frequencies a;i+ and approach each other, as do the corresponding 
families of two- frequency solutions which aX M/C ~ 18.7 coincide at c^b ~ 6.8 
(< uj\- ~ 7.1). At this point they bifurcate with each other, and compared to the 
scenario in figure El (c) an additional simultaneous collision at /3 = +1 appears for 
both solutions, in the regime where for larger M jC the black solution was unstable 
and the grey stable. For smaller M jC the solutions split into a low-frequency and a 
high-frequency part separated by a forbidden gap in cUb as illustrated by figure IHl Thus, 
only the high-frequency part is now connected to the stationary solution. 

Decreasing further M jC the high-frequency part shrinks as ci;i+ and approach 
each other, and furthermore at M jC ~ 18.2 (figure d (a), (c), (e)) the continuation of 
the grey solution from uj^ = changes direction to occur towards larger o^b. Thus, for 
18 < M /C < 18.2 the high-frequency branch consists of one single family which through 
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Figure 7. Continuation of the high-frequency branches of figure|B|(a)-(c) for N = 18.2 
(left figures) and N = 18.01 (right figures). At A/" = 18.2 there is still a tiny remainder 
of the grey part of previous figures which continues approximately 10~^ in the direction 
of negative Wb from the linear value cjb = ~ 7.186 6898. 
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a monotonous continuation versus c^b connects the linear eigenmode at uj\^ with that 
at uji^ (figure m (b), (d), (f)). The solutions are linearly stable close to c^b = but 
unstable for intermediate frequencies (figure [3 (f)) (note also the 'inverse N-shape' of 
if(co'b) in figure[7|(d) consistent with the negative Krein signature of u}\_ and the positive 
of u}\+). At the bifurcation point M /C = 18, ui- = uji+ and the loop of surrounding two- 
frequency solutions shrinks to zero. Consequently, for M /C < 18 where the stationary 
solution is unstable, no nearby two-frequency solutions exist, and the scenario is that of 
a type II HH bifurcation |2S1 (cf. figure 3 in [2H] and figure 2(a) in j2S], [201 )• However, 
the low-frequency branches in figure IHl exist also for M /C < 18, but at a considerable 
distance from the stationary solution (e.g. at A/" = 18 and C = 1 the Hamiltonian (jH]) 
for the stationary solution is if = —27, while its maximum value on the disconnected 
branch of two-frequency solutions is H ^ —39.4). We illustrate the consequences for 
the dynamics of the unstable stationary solution in section ^ 

3.2.3. Smaller M jC . Continuing the low-frequency branches of figure IHl towards 
smaller M /C , their features remain qualitatively the same. The example for M jC = 15 
in figure |H1 shows that they are distinct from those originating from the low-amplitude 
bifurcation (cf. figure EJ. Moreover, they remain at a large distance from the unstable 
stationary solution (jlj) also when J\f/C decreases (see e.g. the large difference in H in 
figure IHl (b)), and their maximum value of c^b (i.e. where the black and grey branches 
meet) continues to decrease (e.g. at A/" = 10 and C = 1 they meet at Uh ~ 2.33). 

To illustrate the different nature of the solutions in figures |3] and |H1 for small cUb, 
we plot in figure IHl the resulting dynamics for one half-period Tb/2 (remembering that 
(piit + T]^/2) = 4>2(t), 4>3(t -|- Tb/2) = (psit) yields the complete dynamics) at cUb = 0.58. 
As C(Jb decreases the fiat parts of the curves will fiatten and stretch out in time, so 
that as tUb — * the previously described asymptotic solutions are approached. (Since 
these solutions are strongly unstable, their accurate numerical determination even at 
cUb = 0.58 and for only half a period requires working in quadruple precision.) 

4. Instability-generated dynamics 

Let us now discuss the dynamics resulting from initial conditions taken as a slightly 
perturbed stationary solution ^ in various parts of its Krein instability regime, with 
three examples shown in figure E3 Close to the low-amplitude instability threshold 
(figure Uni (a)), the dynamics remains trapped forever in a state of bounded oscillations 
around the initial solution. In particular, the amplitude of the initially unexcited site 
{n = 3) always remains small compared to the other two sites. As M /C increases the 
amplitude of these oscillations increases, and above a certain value M jC ~ 10.6 the self- 
trapping is destroyed and the dynamics of all three sites becomes qualitatively the same. 
Slightly above this value the dynamics is still initially trapped (figure [TUKb)), but finally 
the trapping is destroyed and the amplitude of the third site increases drastically. The 
length of the initial transient approaches infinity as M jC 10.6, but decreases for 
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Figure 8. Continuation of the low- frequency branclies of figure IHl for Af — The 
eigenvalues not shown in (c) for small o^b increase monotonously for both solutions 
(/3 ~ 5 • 10^ at Wb = 0.8). As before, horizontal lines in (a), (b) and dashed line in (d) 
show corresponding quantities for the (unstable) stationary solution (0J (C = 1). 



increasing M jC so that close to the high-amphtude instabihty threshold M jC = 18 no 
transient is observed (figure ITIll (c)). In all cases, the dynamics in the untrapped regime 
can be characterized as 'intermittent population inversion', which was observed also for 
certain initial conditions of the (nonequivalent) trimer studied in This is most 

clearly seen in figure EH (c). Here, one can during certain time intervals unambiguously 
identify two sites of large amphtude and a single small-amplitude site, which changes 
from initially being n = 3 to n = 2, and then, during a rather long time, to = 1. After 
that follows a regime of more complicated oscillations, which then again is followed by 
a regime with a single small- amplitude site changing from n = 3 to n = 1, etc. 

To understand more clearly the origin of this self-trapping transition, we analyze the 
dynamics close to the unstable stationary solution in terms of Poincare sections. They 
can be introduced in many different ways; here we apply similar ideas as in jn2] making 
use of the transformation into action-angle variables {P„, 6'„} defined by ipn = y/Pn^~^^" 
(a slightly different approach was used in It is then convenient to replace one 

of the action variables, which we here choose as Pi, with the conserved quantity Af. 
Then, the angle- variables conjugated to the set of generalized momenta {A/", P2, P3} are 
|n2] {01,02 — 01,03 — 0i}, so that 01 describing an overall phase becomes an ignorable 
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Figure 9. Dynamics for one half-period of the four different solutions in figures^ 
and Oat Af — 15, — 0.58 and C — \. (a) and (c) [(b) and (d)] correspond to the 
black [grey] solutions in figures ^ and |S1 respectively. 



coordinate. Thus, the essential dynamics takes place in a four-dimensional space where 
the surface of constant energy H is three-dimensional, so that a proper Poincare section 
through it becomes two-dimensional. Consequently, although chaotic trajectories may 
fill a large portion of the available phase space [2S1 112] , Arnold diffusion is prohibited 
|3H EH] since the presence of any regular KAM tori will disconnect the phase space. 
(Arnold diffusion does however appear for the four-site DNLS model l5Uj .) 

As a particular choice of Poincare section giving a clear illustration of the self- 
trapping transition, we plot in figure ITD = l^/'sp versus ^3 — 9i at each time instant 
when 92 — 9i = 7c and ^{02 — di) < 0. Since the stationary solution has zero amplitude at 
77, = 3 its phase 63 is undefined, and thus it is represented by a horizontal line iV'sP = 0. 
The scenario in the self-trapped regime is illustrated by figure ITTl (a), corresponding to 
the dynamics in figure ITUl (a). Here the elliptic fixed point at ^3 — 6*1 = vr represents the 
stable two-frequency solution of section IXD with the same value of H as the stationary 
solution. As in figure El (d) this solution belongs to the black branch. This elliptic fixed 
point appears at I'lpsl'^ = at the type I HH bifurcation point M /C ~ 9.077, and moves 
vertically in the direction of increasing for increasing M /C (cf. figure El (b)). It is 
surrounded by two different kinds of regular periodic or quasiperiodic orbits, where the 
latter constitute KAM tori corresponding generically to quasiperiodic i/iree-frequency 
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solutions in the original DNLS dynamics. With a pendulum analogy, these orbits can 
be classified as 'rotating' and 'vibrating', respectively, where the former extend for all 
values of 6*3 — 61 while the latter only exist in a bounded region close to ^3 — 6*1 = vr. 
Then, the unstable stationary solution becomes the separatrix between these different 
kinds of solutions. Although the separatrix is chaotic (which can be seen from a careful 
look at figure ^2 (a)) the chaos is confined between KAM-tori, and in particular the 
existence of confining 'rotating' tori makes it impossible for iV'sP to exceed some upper 
limit value (iV'sP ~ 1 at 6^3 — 6'i = tt in figure ITT] fa)). Thus, self-trapping results. 

Below M jC ~ 9.077, where the stationary solution is stable and the elliptic 
fixed point with nonzero |V^3p not yet born, all surrounding KAM tori are of the 
'rotating' kind. As M jC is increased and the elliptic fixed point moves upwards towards 
larger iV'sP? more and more of the 'rotating' KAM tori get destroyed, and finally at 
M jC ~ 10.6 the last 'rotating' KAM torus breaks up, and the self-trapping is destroyed. 
The Poincare plot for M jC = 10.7, i.e. just above the transition point, is shown in figure 
ITT] (b) (the corresponding dynamics is qualitatively similar to figure HU] (b) but with a 
longer transient t ~ 17000). It can be seen, that although the dynamics finally spreads 
to a large part of the available phase space, the darker parts signifies regions where it 
will be almost trapped for long times. This should be expected close to the transition 
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— TT, ^(6*2 ^ 6*1) < for the dynamics close to 
unstable stationary solutions Q with C = 1 and (a) TV = 10 (i? = 5), (b) TV = 10.7 
[H = 3.4775), and (c), (d) TV = 17 (ff = -21.25), respectively. The stationary solution 
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0. The elliptic fixed point (filled circle) 



represents continuations of the stable branch of two- frequency solutions in figures EE 
The window of regular orbits in (d) is marked by a short vertical black line in (c). 



point, as destroyed KAM tori should transform into cantori by a transition by breaking 
of analyticity jSI]. Close to the transition the cantori should have almost full measure, 
and can thus trap the dynamics for very long times. (The white regions in the upper 
part around ^3 — = mod 27r correspond to dynamics trapped around the equivalent 
two- frequency solution with main amplitudes at sites n = 2 and n = 3.) 

At the self-trapping transition M /C ~ 10.6, the solution family represented by the 
elliptic fixed point in figure ^2 changes from the black (larger H) to the grey (smaller 
H) branch (cf. section IB. 1.31 and figure El (b)). Increasing further M jC, this fixed point 
continues to move towards larger |'?/'3p (cf. figure IH (a)), and the surrounding island 
of KAM tori shrinks so that the unstable stationary solution invades almost all the 
available phase space (figure [TT] (c), (d), corresponding to the dynamics of figure ITIH 
(c)). Here, there are no visible traces of destroyed rotating KAM tori (if cantori exist 
they should have very small measure), so the dynamics shows no transient trapping. 
This behaviour persists until the stationary solution regains stability through the type 
II HH bifurcation at TV/ C = 18. Before, at TV/ C ~ 17.2, the grey two-frequency solution 
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becomes unstable through a period-do ubhng type bifurcation (cf. colhsion at /3 = — 1 
in figure E](c)), so that the eUiptic fixed point in figure ITT] fd) becomes hyperbohc and 
sphts into two new elhptic points surrounded by tiny KAM tori. 

5. Conclusions 

In conclusion, we have investigated the Hamiltonian Hopf bifurcations in the three-site 
DNLS model, resulting in a regime of oscillatory instability for the stationary solution 
with two sites of anti- phased oscillations and the third site of zero amplitude ('single 
depleted well' Using numerical continuation techniques to calculate the (generally 

quasiperiodic) two-frequency solutions involved in the two bifurcations, we found them 
to be of two different types. At the low-amplitude instability threshold the bifurcation is 
of 'type r, with stable two- frequency solutions close to the unstable stationary solution 
near the bifurcation point. These two- frequency solutions are themselves surrounded 
by KAM tori, representing generally quasiperiodic three-frequency solutions of the full 
DNLS model. At the high-amplitude instability boundary, the HH bifurcation is of 
'type ir, and the stationary solution has no surrounding two- frequency solutions on 
the unstable side of the bifurcation. This reflects itself as a self-trapping transition 
in the instability-induced dynamics of the stationary solution, so that in the low- 
amplitude instability regime the dynamics remains trapped close to the initial state 
with small amplitude on the third site, while in the high-amplitude instability regime 
an 'intermittent population inversion' dynamics is observed with the small-amplitude 
oscillation moving chaotically between the sites. This self-trapping transition, which we 
believe has not been described earlier in the literature, was found to correspond to the 
destruction of phase-space dividing KAM tori. 

As discussed in the introduction, our original motivation was to obtain an 
increased understanding for the dynamics resulting from oscillatory instabilities of 
certain multibreather conflgurations in inflnite lattices, which have been discovered 
recently in many contexts. Based on the results obtained here, we will address these 
issues in a forthcoming publication. However, the three-site DNLS model is of interest 
in itself, in particular in view of the recent experimental progress in the flelds of coupled 
optical waveguides and Bose-Einstein condensates. It seems likely, that the existence 
of stable families of two-frequency solutions, corresponding to intensities periodically 
oscillating around their mean values, as well as the existence of a new type of self- 
trapping transition in the unstable dynamics of the stationary solution, should be 
experimentally verifiable within these contexts. 
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